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ABSTRACT 

This  study  deals  with  the  evaluation  of  monthly  reference  evapotranspiration  (ET0)  estimation  using  FAO-24 
Modified  Penman  (MP)  method  by  comparing  its  performance  with  that  of  FAO-56  Penman-Monteith  (PM)  method, 
developing  relationship  between  MP  and  PM  methods  and  recalibrating  the  method  with  respect  to  PM  method.  Tirupati 
region  of  Andhra  Pradesh,  India  is  selected  as  the  study  area  and  its  meteorological  data  was  collected  from  the  India 
Meteorological  Department,  Pune.  During  comparison,  it  is  observed  that  the  percentage  deviations  of  ET0  values 
estimated  by  MP  method  with  reference  to  PM  method  are  significant.  Even  after  developing  inter-relationships  between 
PM  and  MP  methods,  the  performance  is  not  satisfactory.  Therefore,  it  is  recalibrated  with  respect  to  PM  method. 
The  recalibrated  MP  method  yielded  the  least  RMSE  and  high  R2  &  EC  values  resulting  in  ET0  comparable  with  that  from 
PM  method.  The  MP  method  performed  well  in  terms  of  performance  evaluation  criteria.  The  slope  and  intercept 
respectively  close  to  one  and  zero  also  indicate  an  improved  performance  of  the  method  with  recalibrated  coefficients. 
Therefore,  recalibrated  MP  method  may  be  recommended  for  reasonable  ET0  estimation  in  the  study  area. 
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INTRODUCTION 

Evapotranspiration  (ET)  estimation  is  vital  in  irrigation  system  design,  crop  yield  simulation  and  water  resources 
planning  and  management.  Field  measurement  of  evapotranspiration  is  rarely  available  and  actual  crop  evapotranspiration 
(ETC)  is  usually  calculated  from  reference  evapotranspiration  (ET0)  using  the  crop  factor  method,  which  consists  of 
multiplying  ET0  with  crop  coefficients  (Kc)  to  obtain  ETC  (i.e.,  ETC  =  ET0  x  Kc).  Several  reports  on  the  estimation  of  Kc  are 
available.  Doorenbos  and  Pruitt  (1977)  m,  Doorenbos  and  Kassam  (1979)  m,  Jensen  et  al.  (1990)  [5]  and  Allen  et  al.  (1998) 
[1]  have  reported  crop  coefficients  for  many  crops.  These  values  are  commonly  used  in  places  where  the  local  data  is  not 
available. 

It  is  desirable  to  have  a  method  that  estimates  reasonably  the  ET0.  The  amount  of  ET0  is  calculated  by 
meteorological  data  based  methods.  However,  choosing  the  best  method  for  ET0  estimation  from  the  described  methods  is 
difficult.  Therefore,  the  International  Commission  for  Irrigation  and  Drainage  and  the  Food  and  Agriculture  Organization 
of  the  United  Nations  (FAO)  Expert  Consultation  on  revision  of  FAO  methodologies  for  crop  water  requirements  have 
recommended  the  FAO-56  Penman-Monteith  (PM)  equation,  which  was  presented  in  Allen  et  al.  (1998)[1],  to  be  used  as  the 
standard  method  to  estimate  ET0.  Many  researchers  thereafter  took  it  as  a  standard  to  modify  other  methods  that  required 
less  input  data. 
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Allen  et  al.  (1998)  [11  presented  that  the  Penman  method  may  require  local  calibration  of  the  wind  function  to 
achieve  satisfactory  results.  Kotsopoulos  and  Babajimoponlos  (1997)  [6]  derived  mathematical  expressions  that  describe  the 
parameters  used  for  the  calculation  of  the  Penman  reference  evapotranspiration  through  a  nonlinear  regression  procedure. 
Comparisons  of  these  expressions  to  those  found  in  the  literature  showed  more  reliable  results.  Hussein  (1999) [4]  compared 
Penman  combination  equations  with  measured  ET  data  from  a  clipped  grass  lawn.  The  PM  method  gave  the  best 
agreement  with  measured  data,  followed  by  the  FAO  Penman  and  Penman-Watts-Hancock  methods.  The  1963  Penman 
method  was  the  least  satisfactory  among  the  methods  tested.  However,  its  performance  improved  significantly  when  the 
saturation  vapour  pressure  was  calculated  as  the  mean  value  at  maximum  and  minimum  air  temperatures.  Shah  and  Edling 
(2000)  191  evaluated  PM,  FAO-Penman  and  1963  Penman  combination  models  for  their  capabilities  to  predict  rice  ET  using 
daily  weather  data.  The  PM  (daily)  method  had  a  R2  of  63.7%  as  compared  with  62.9,  60.0  and  61.7%  for  PM  (hourly), 
FAO  Penman  and  1963  Penman  methods  respectively. 

The  present  study  reports  the  performance  evaluation  of  MP  method  based  on  its  accuracy  of  estimation  and 
development  of  inter-relationship  between  MP  and  PM  methods  and  also  recalibrated  with  reference  to  PM  method  for 
Tirupati  region  of  Andhra  Pradesh,  India. 

MATERIALS  AND  METHODS 

Tirupati  region,  located  in  Chittoor  district  of  Andhra  Pradesh,  India,  with  global  coordinates  of  13°  05' N  latitude 
and  79°  05'E  longitudes,  has  been  chosen  as  the  study  area.  The  meteorological  data  of  the  study  area  for  the  period 
1992-2001  was  collected  from  IMD,  Pune.  Data  from  1992  to  1998  is  used  for  the  purpose  of  training  the  model  and  that 
of  1999  to  2001  for  testing  the  model.  The  details  of  the  methods  selected  for  the  present  study  are  presented  in  Table  1. 


Table  1:  Details  of  Reference  Evapotranspiration  Estimation  Methods 


Method 

FAO-56 

Penman- 

Monteith 

(PM) 

Method 

Allen  et  al., 
(1998)[1] 

0.408A  (Rn  -  G ' )  +  r       9°°      u2  (es  -  ea ) 

j^rp   _  mean 

A  +y  (l  +  0.34«2) 

T  T 

1  max,  1  min, 

RHmax, 
RHmin,  U2, 

n 

FAO-24 

Modified 

Penman 

(MP) 

Method 

Doorenbos 
and  Pruitt 

(1977) [3] 

ET0  =  C  _^_R  +^L_(0.27)(i.o  +  0.01U2)(e,-e.)] 
|_A  +  y        A  +  y 

where 

C  =  0.68  +  0.0028  (RHmax)  +  0.018  (Rs)  -  0.068  (ud) 
+  0.013  (ud  /  un)  +  0.0097  (ud)(ud/un) 
+  0.000043  (RHmax)  (Rs)  (ud) 

T  T 

1  max,  1  min, 

RHmax, 
RHmin,  n 

u2,  ud/un 

PERFORMANCE  EVALUATION  CRITERIA 

The  performance  evaluation  criteria  used  in  the  present  study  are  the  coefficient  of  determination  (R2),  the  root 
mean  square  error  (RMSE),  systematic  RMSE,  unsystematic  RMSE  and  the  efficiency  coefficient  (EC). 

Coefficient  of  Determination  (R2) 

It  is  the  square  of  the  correlation  coefficient  (R)  and  the  correlation  coefficient  is  expressed  as 
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R  = 


X      (°i  ~  °)(Pi  ~  P) 

 i  =  l  

t  (o.-oft  (p.-p-y 


Where  O  and  P  are  observed  and  estimated  values,  ^  and  P  are  the  means  of  observed  and  estimated  values  and 
n  is  the  number  of  observations.  It  measures  the  degree  of  association  between  the  observed  and  estimated  values  and 
indicates  the  relative  assessment  of  the  model  performance  in  dimensionless  measure. 

Root  Mean  Square  Error  (RMSE) 

It  yields  the  residual  error  in  terms  of  the  mean  square  error  and  is  expressed  as  (Yu  et  al.,  1994) [101 


RMSE  = 


n 


Systematic  RMSE  (RMSES) 

It  measures  the  room  available  for  local  adjustment.  It  is  expressed  as 


RMSE  = 


Z  (A-«,)2 


;=-] 


n 


Where  1 '  ' ,  a  and  b  are  the  liner  regression  coefficients 

Unsystematic  RMSE  (RMSEU) 

It  shows  the  noise  level  in  the  model  and  is  a  measure  of  scatter  about  the  regression  line  and  potential  accuracy. 
It  is  expressed  as 

t  (Pi-Pi)2 

J=l  

n 

Efficiency  Coefficient  (EC) 

It  is  used  to  assess  the  performance  of  different  models  (Nash  and  Sutcliffe,  1970)  .  It  is  a  better  choice  than 
RMSE  statistic  when  the  calibration  and  verification  periods  have  different  lengths  (Liang  et  al.,  1994)  .  It  measures 
directly  the  ability  of  the  model  to  reproduce  the  observed  values  and  is  expressed  as 


RMSEU  = 
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EC  =  1--^- 


A  value  of  EC  of  90%  generally  indicates  a  very  satisfactory  model  performance  while  a  value  in  the  range 
80-90%,  a  fairly  good  model.  Values  of  EC  in  the  range  60-80%  would  indicate  an  unsatisfactory  model  fit. 

RESULTS  AND  DISCUSSIONS 

Comparing,  developing  inter-relationship  and  recalibrating  the  MP  method  with  respect  to  PM  method. 

Comparison  between  MP  Method  and  PM  Method 

The  percentage  deviations  of  ET0  values  estimated  by  MP  method  with  reference  to  PM  method  are  11.1  to  36.2. 
The  deviations  are  significant.  Fig.l  showing  the  comparison  of  ET0  estimates  with  those  of  PM  ET0  also  exhibit  similar 
observations 
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Figure  1:  Comparison  of  Monthly  Average  ET0  Values  Estimated  by  MP  Method  with  Reference  to  PM  Method 
Development  of  Inter-Relationship  between  MP  Method  and  PM  Method 

The  ET0  values  estimated  using  MP  method  plotted  against  PM  ET0  are  shown  in  Figure  2.  The  performance 
indicators  of  the  relationship  developed  between  PM  and  MP  methods  are  presented  in  Table  2.  Reasonably  high  R2  values 
indicate  a  strong  correlation  between  PM  ET0  and  MP  ET0  values.  This  may  be  due  to  the  fact  that  both  the  methods  are 
physically  based  (combination)  and  require  similar  climatic  parameters  for  ET0  estimation.  But  large  deviations  of  the 
slope  from  one  and  intercept  from  zero  indicate  the  unsatisfactory  performance  of  relationships  between  PM  and  MP 
method.  And  also,  the  considerable  magnitude  of  unsystematic  RMSE  indicates  the  requirement  of  further  dampening  of 
noise  level  in  the  method  represent  the  scope  for  recalibration.  Therefore,  it  may  be  tried  to  improve  this  method's 
performance  by  suitably  recalibrating  it  against  PM  method  using  the  observed  climatic  data. 
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Figure  2:  Scatter  Plot  of  Monthly  Average  ET0  Values  Estimated  by  MP  Method  against  PM  Method 
Impact  Factor  (JCC):  5.7179  Index  Copernicus  Value  (ICV):  3.0 
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Method 

Relationship 

R 

RMSE  (nun) 

RMSEs(mm) 

RMSEu  (nun) 

EC  (%) 

MP 

PM  =  0.8405MP  -  0.2833 

0.9536 

0.29 

0.06 

0.29 

95.36 

Recalibration  of  MP  Method  with  Reference  to  PM  Method 

It  has  been  emphasized  in  the  above  section  that  MP  method  selected  for  the  present  study  has  not  performed 
satisfactorily  in  the  regional  ET0  estimation.  The  relationship  developed  between  PM  method  and  this  method  to  estimate 
ET0  showed  an  unsatisfactory  performance  in  terms  of  RMSE.  Therefore,  it  is  recalibrated  with  respect  to  PM  method. 

The  recalibrated  equation  derived  for  the  study  area  is  presented  along  with  original  equation  in  Table  3. 
The  performance  indicators  of  recalibrated  equation  in  the  estimation  of  ET0  for  both  training  and  testing  periods  are  given 
in  Table  4.  The  scatter  and  comparison  plots  of  ET0  values  estimated  by  this  method  with  reference  to  PM  ET0  during  the 
testing  period  are  shown  in  Fig. 3  and  Fig.4  respectively. 

It  may  be  observed  from  Table  4  that  the  MP  method  yielded  the  least  RMSE  and  high  R2  &  EC  values  resulting 
in  ET0  comparable  with  that  from  PM  method.  The  regression  coefficients  of  scatter  plots  close  to  one  and  zero  also 
indicate  considerable  dampening  of  noise  level  in  the  equation  with  recalibrated  coefficients  by  decreasing  the  scatter 
about  the  regression  line. 

From  the  above  discussion,  it  may  be  concluded  that  the  MP  method  with  recalibrated  coefficients  may  be  used 
for  reasonable  ET0  estimation  in  the  study  area. 

Table  3:  Recalibrated  MP  Equation 


Method 


FAO-24 

Modified 

Penman 

(MP) 

Method 


Original  Equation 


ET0  = 
C 

A 


R„  +- 


A  +  y 


A  +  y 
Where 

C  =  0.68  +  0.0028  (RHmax) 
+  0.018  (Rs) 

-0.068  (ud)  + 0.013  (ud/un) 

+  0.0097  (ud)(ud/u„) 

+  0.000043  (RHmax)  (Rs)  (ud) 


(0.27)(1.0  +  0.01U,)(es-eJ 


Recalibrated  Equation 


ET0  = 
C 

A 


R„+- 


(0.27)(1.0+0.01U2)(es-eJ 


A+y  "  A+y 
Where 

C  =  0.96 -0.001 19  (RHmax) 
-0.0108  (Rs) 

-0.045  (ud)  +  0.013  (ud/u„) 

+  0.0097  (ud)(ud/un) 

+  0.000032  (RHmax)  (Rs)  (ud) 


Table  4:  Performance  Indices  of  Recalibrated  MP  Method 


Method 

Slope  of  the 
Scatter  Plots 

Intercept  of  the 
Scatter  Plots 

R2 

RMSE  (nun) 

EC  (%) 

Training 
Period 

Testing 
Period 

Training 
Period 

Testing 
Period 

Training 
Period 

Testing 
Period 

Training 
Period 

Testing 
Period 

Training 
Period 

Testing 
Period 

MP 

1.0084 

1.0077 

-0.0420 

-0.0490 

0.9976 

0.9970 

0.07 

0.07 

99.76 

99.70 
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Figure  3:  Scatter  Plots  of  Monthly  Average  ET0  Values  Estimated  by  Recalibrated 
MP  Method  against  PM  ET0  Values  during  Testing  Period 
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Figure  4:  Comparison  of  Monthly  Average  ET0  Values  Estimated  by  MP  Method  with 
Reference  to  PM  Method  during  Testing  Period 

CONCLUSIONS 

The  percentage  deviations  of  ET0  values  estimated  by  MP  method  with  reference  to  PM  method  are  significant. 
The  MP  method  has  not  performed  well  even  after  developing  inter-relationship  with  PM  method.  The  MP  method 
improved  its  performance  significantly  on  recalibration.  The  recalibrated  MP  method  performed  equally  well  with  PM 
method.  Therefore,  recalibrated  MP  method  may  be  recommended  for  reasonable  ET0  estimation  in  the  study  area. 
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